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PREFACE 

This  report  documents  an  investigation  at  Rand  of  the  laminar 
boundary  layer  over  heated  bodies  undergoing  either  acceleration  or 
deceleration  in  water.  The  research  is  part  of  a larger  Rand  re- 
search project,  sponsored  by  the  Defense  Advanced  Research  Projects 
Agency,  which  is  engaged  in  developing  hydrodynamic  design  techniques 
that  exploit  various  methods  of  boundary- layer  control.  Tlie  report 
presents  a method  of  predicting  laminar  boundary-layer  characteristics 
when  heating,  pressure  gradients,  and  degrees  of  unsteadiness  are 
acting  in  concert.  Using  these  results  and  boundary-layer  stability 
theory,  the  enhanced  boundary-layer  stability  produced  by  an  acceler- 
ating flow  is  illustrated. 

The  report  should  be  useful  to  hydrodynamiclsts,  designers  of 
submersibles,  and  others  engaged  in  fluid  mechanics  research.  Other 
related  Rand  publications  include: 

R-1752-ARPA/ONR,  Low-Speed  Boundary -Layer  Transition  Workshop, 
June  1975. 

R-1789-ARPA,  Controlling  the  Separation  of  Laminar  Boundary 
Layers  in  Water:  Heating  and  Suction,  September  1975. 

R-1863-ARPA,  The  Effects  of  Wall  Temperature  and  Suction  on 
Laminar  Boundary -Layer  Stability,  April  1976. 

9 

R-1898-ARPA,  "e  Stability  Theory  and  Boundary-Layer  Transi- 
tion, February  1977. 

R-1966-ARPA,  The  Buoyancy  and  Variable  Viscosity  Effects  on  a 
Water  Laminar  Boundary  Layer  Along  a Heated  Longitudinal 
Horizontal  Cylinder,  February  1977 . 

R-2111-ARPA,  Entry  Flow  in  a Heated  Tube,  June  1977. 

R-2165-ARPA,  Approximate  Methods  for  Calculating  the  Properties 
of  Laminar  Boundary  Layers  in  Water  (to  be  published). 
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SUMMARY 


The  laminar  boundary  layer  over  heated  bodies  undergoing  acceler- 
ation or  deceleration  In  water  has  been  investigated  to  determine  the 
combined  influence  of  heating,  pressure  gradient,  and  unsteadiness  on 
boundary- layer  stability.  For  simplicity,  the  boundary- layer  flow 
over  accelerating  wedges  or  cones  has  been  analyzed  for  values  of  the 
acceleration  parameters,  ■ (x/l)|) (dU^/dt) , which  are  small  but  not 

negligible.  It  is  shown  that  the  boundary- layer  stability  analysis, 
in  this  range  of  e^,  requires  solution  of  the  unsteady  boundary- layer 
equations,  but  that  stability  characteristics  may  be  determined  using 
the  classic  Orr-Sommerf eld  approach. 

Velocity  profiles  have  been  calculated  by  a perturbation  expan- 
sion in  about  the  quasi-steady  flow,  and  the  stability  of  the  result- 
ing profiles  has  been  calculated  by  an  accurate  and  convenient  modifi- 
cation of  the  Dunn-Lln  theory,  specialized  to  constant  density,  variable 
property  flow.  The  results  confirm  that  the  acceleration  is  stabilizing, 
while  deceleration  is  destabilizing,  and  show  that  the  stability  enhance- 
ment of  acceleration  is  less  dramatic  when  both  positive  velocity  gra- 
dient and  heating  are  present.  A value  of  of  only  .1  produces  a four- 
fold increase  in  the  critical  Reynolds  number  for  a flat  plate  and  a 
two-fold  increase  for  a highly  favorable  pressure  gradient.  Heating 
reduces  these  changes  by  as  much  as  a factor  of  two.  The  probable  ex- 
planation of  the  effect  Is  discussed. 
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I . INTRODUCTION 


Unsteady  laminar  boundary  layers  are  found  in  a wide  variety  of 
practical  circumstances,  but  the  quantitative  evaluation  of  the  un- 
steady effects  is  not  yet  widely  considered.  The  effects  of  unsteadi- 
ness are  usually  assessed  by  employing  the  quasi-steady  approximation. 
However,  for  a slightly  more  vigorous  unsteadiness,  the  calculation 
procedure  has  been  to  consider  small  departures  from  the  quasi-steady 
state.  This  theory  is  well  developed  and  is  discussed  in  Refs.  1 
through  3 for  arbitrary  magnitudes  of  slowly  varying  unsteadiness.  The 
accuracy  of  this  theory  is  presently  estimated  from  order  of  magnitude 
estimates,  or  from  the  calculation  of  further  terms  in  the  solution 
sequence.  However,  a recent  theoretical  study  of  fluctuating  free- 
stream  perturbations  to  flat-plate  flow  via  matched  asymptotic  methods 
suggests  strongly  tliat  the  quasi-steady  approach  is  uniformly  valid 
when  the  frequency  or  flow  length  is  small  in  comparison  to  a reduced 
frequency  or  a body  length,  respectively  (Ref.  4).  Nevertheless,  it 
appears  that  an  assessment  of  the  accuracy  of  the  technique  must  await 
ultimate  comparison  with  results  from  the  newly  fledged  numerical 
schemes  for  unsteady  flow.  Until  this  is  done,  the  quasi-steady  approx- 
imation sequence  used  here  must  be  considered  the  practical  method  of 
choice  for  realistic  mild  unsteadiness  (Ref.  5). 

This  report  discusses  the  effects  of  unsteadiness  on  laminar 
boundary- layer  flows  and  their  stability  as  governed  by  the  Orr- 

"k 

Sommerfold  equation.  It  presents  tlie  required  conditions  for  the 
validity  of  the  quasi-steady  approximation  to  be  a valid  approximation 
of  the  governing  flow  equations.  It  is  well  known  that  a critical 
parameter  for  boundary- layer  flow  is  the  ratio  of  diffusion  time  to 
flow  time;  the  criterion  for  nearly  quasi-steady  flow  Is  that  this 
parami'ter  should  he  small.  It  is  also  shown  that  t lu'  criteria  for 

* 

U.  Von  Kercxi'k  of  Naval  Ship  Researi’h  ami  Development  th'iiter 
(NSRDC)  discussed  the  effect  on  transition  In  a presentation  at  the 
Rand  Workshop  on  Low  Speed  Koundary-I.ayer  Transition.  Santa  Mc'nica, 

('.a  Ilf.,  September  1976. 
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valldlty  of  the  Orr-Sommerfeld  equation  are  based  on  modifications 
of  this  ratio,  and  are  not  nearly  as  stringent.  Examples  of  heated 
wedge  flows  in  water  are  presented  and  discussed  that  show  the  pro- 
found effect  of  even  slowly  varying  unsteadiness  on  both  laminar 
boundary-layer  flow  and  its  stability. 


I 

i 

I 


i 


) 

( 
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II.  ANALYSIS 


Vfhen  Che  flow  external  Co  Che  boundary  layer  varies  with  time, 
variations  in  boundary-layer  momentum  are  produced,  diffused  to  the 
wall,  and  dissipated.  It  is  through  this  mechanism  that  Che  boundary- 
layer  flow  accommodates  variation  in  free-stream  velocity.  If  the 
boundary  layer  accommodates  instantaneously,  the  flow  may  be  considered 
quasi-steady.  If,  on  Che  other  hand,  the  boundary  layer  accommodates 
slowly,  Che  flow  is  unsteady  aud  the  rate  of  accommodation  is  important. 
\ measure  of  how  rapidly  the  boundary  layer  accommodates  is  the  ratio 
of  diffusion  time  to  flow  time.  In  Che  following  discussion  this  is 
the  critical  parameter,  and  all  comparisons  of  unsteadiness  in  boundary- 
layer  flow  to  unsteadiness  in  boundary-layer  stability  are  based  on  its 
magnitude. 

The  analysis  will  first  focus  on  bound ary- layer  flow  in  that  the 
quasi-steady  flow  parameters  will  be  derived  and  numerical  examples 
will  be  presented,  and  next  on  boundary-layer  stability  in  that  a cri- 
terion for  nearly  quasi-steady  flow  will  be  developed  and  numerical 
examples  of  Che  effect  of  unsteadiness  on  critical  Reynolds  numbers 
will  be  shown.  The  specific  flows  Chat  are  used  as  examples  are  two- 
dimensional  Falkner-Skan  flows  with  wall  heating,  variable  viscosity, 
and  Prandtl  number  corresponding  to  water. 


BOUNDARY-LAYER  FLOW 

The  unsteady  equations  for  a two-dimensional  laminar  boundary- 
layer  flow  at  an  incompressible  fluid  with  heat  transfer  and  variable 
viscosity  are 


3u  3u 

9t  3x 


+ 


_l£_  + _iL 

p3x  p3y 


I? 
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(1) 


3T  ^ 3T  ^ 3T  '^e  3^T 

3t  ^ 3x  3y  Pr  ^^2  ’ 


(2) 


where  Pr  is  based  on  freestream  properties. 
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To  assess  tlio  etfei-t  of  unsteadiness  on  Inuindarv- 1 ayer  flow,  the 
relative  order  of  magnitude  ol  tlie  first  term  on  tlie  lelt-liand  side 
of  tile  momentum  equation  and  ^lie  viscous  term  on  llte  riglit  side  must 
be  determined.  The  order  of  magnitude  of  the  viscous  term  is 


3y 


wliere  Uj  is  the  external  velocity  with  an  arbitrary  time  dependence. 


and  6 is  the  momentum  layer  tlilckness. 

Tlie  order  of  magnitude  of  the  unste.-id iness  is  estimated  by  tlie 
variation  in  time  of  the  freest  ream  velocity. 


Ti  * ^at”)  • 


Then,  the  ratio  of  magnitude  of  the  two  terms  is 


6 ^ unsteady  effect 

V 3tUj  viscous  effect 


A more  physical  interpretation  consistent  witli  that  described  in 
the  Introduction  is  obtained  by  noting  tliat 


3U 


1 

3tir, 


where  t^.  Is  a characteristic  flow  time,  i.e.,  the  time  it  takes  tlie 


external  flow  to  change  from  one  state  to  another.  Further,  v/6  is 

,2 


a diffusion  velocity  and  6 /v  is  the  time  it  takes  the  momentum  to 
diffuse  to  the  wall.  Therefore,  the  ratio  of  the  magnitude  of  the 
unsteady  effect  to  the  viscous  effect  is  related  to  the  ratio  of  char- 
acteristic times. 
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V 3tU. 


If  the  bt)uniliiry- layer  thickness  Is  related  to  the  Incal  Reynolds 
_c 

number  [6  ~ x(Re)  ^|,  the  quasi-steady  parameter  of  Refs.  1 through  3 
Is  obtained: 


.2  3U,  3U, 

o_ ^ X 1 

V 3tU'  ~ 3t  ' “ 1 ■ 


(3) 


When  this  parameter  is  small,  or  1»  the  flow  in  the  boundary 

layer  is  quasi-steady.  A similar  parameter  is  discussed  in  Ref.  4. 

For  flows  that  are  nearly  quasi-steady,  a systematic  series  expansion 
procedure  can  be  developed.  A detailed  discussion  of  this  procedure 
as  it  applies  to  general  unsteady  flow  is  given  in  Ref.  2.  However, 
a special  version  of  nearly  quasi-steady  flow  and  this  technique  are 
illustrated  in  the  example  discussed  below. 

Unsteady  Boundary- Layer  Flow  over  Heated  Bodies 

The  theory  developed  In  Refs.  2 and  3 will  be  extended  to  examine 
the  unsteady  effect  on  a flow  over  heated  bodies  in  water  with  constant 
freestream  and  wall  temperatures.  Although  the  general  theory  will  be 
developed,  this  discussion  is  concerned  with  nearly  quasi-steady  flows 
where  the  time  variation  is  limited  to  slowly  varying  functions  and 
Wedge  flows.  The  attractive  attribute  of  this  theory  is  that  for  nearly 
quasi-steady  flows  the  unsteady  i-otnponent  of  the  freestream  velocity  can 
have  arbitrary  magnitude. 

Following  Ref.  ),  a boundary- 1 ayer  scaling  will  bo  employed  that 
allows  the  quasi-steady  component  to  be  iletermlneil  .is  the  zeri'-order 
solution.  I’he  sc.iling  to  transform  Kqs.  (1)  .and  (2)  is 
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\/^ 

V 2vx 


<p  • ^2vxUj  f(n,x,t)  , 


(4) 


T - T 
e 

T - T ’ 
w e 


where  U^,  the  Inviscid  unsteady  velocity  along  the  surface.  Is  given 
by  - Ug(x)  A(t). 

Note  that  this  form  of  the  time-dependent  inviscid  surface  veloc- 
ity corresponds  to  unsteady  irrotational  flow,  in  the  absence  of  dis- 
placement effects.  V!g(x)  corresponds  to  the  surface  velocity  when  a 
body  translates  with  unit  velocity  into  a fluid  at  rest  and  A(t)  is 
the  actual  time-dependent  velocity  of  translation.  Thus,  Ug(x)  is  to 
be  determined  from  a conventional  potential  flow  solution  corresponding 
to  unit  body  velocity.  The  transformed  equations  for  the  reduced  stream 
function,  f(x,r|.t),  and  temperature,  0(x,r|,t),  are 


= 2x(f  f 

n nx 


2U  X 

*nn'x>  * ^ 1 'nn 


1) 


(5) 


and 


+ ^ f 

nt  ’ 


e \ 

-r^  + ( 1 + M(x))  f0  = 2x(f  0 - 0 f ) 

Pr  \ / n n X n X 


+ — 0 ) + — 0 
y2  4 ^ '^t  ’ 


(6) 


F 


where  M(x)  ratio  of  boundary  layer  to  free- 

stream  viscosity.  N can  be  approximated  by  the  empirical  expression 
for  water  used  in  Ref.  6: 


7 1 ^ -1 

N - y/Ug  • (a  + br  + cr  + dr  + er  ) -^j — , (7) 

e 

where  r » T/491.69®R,  and  Uj.g£/Ug  * 4.339  Ib/hrft/y^.  The  appropriate 
boundary  conditions  are 

t ^ 0 , 

n**’”®,  f^-*"!,  0-^O,  (8) 

T1  = 0,  f = f^  = 0,  0-1. 

The  general  velocity  and  temperature  distributions  are  functions 
of  three  Independent  variables,  l.e., 

f = f(n,x,t)  . (9) 

However,  for  special  cases  where  f is  a slowly  varying  function  of 
time,  one  will  find  it  convenient  to  express  the  velocity  distribution 
as 


f = f(Ti,x;e^)  , 


(10) 


where  now  f depends  on  a sequence  of  general  unsteady  parameters 


e 

n 


ft)  V' 


n = 1,2,3  . . . 


(11) 


The  first  parameter,  is  an  obvious  expansion  parameter  (the  ratio 

of  diffusion  time  to  flow  adjustment  times)  that  is  explicitly  displayed 
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in  the  equations.  For  n > 1,  tile  parameters  are  ratios  of  various 

diffusion  times  to  various  flow  times  and  a result  of  suhstiiutinn  ll'e 

expansion  of  f into  the  governing  equations  and  requiring  a self- 

consistent  set  of  differential  equations. 

When  a body  has  constant  acceleration,  the  parameters  c for 

n 

n • 1 are  precisely  zero.  For  flows  that  have  nearly  lonstant  accel- 
eration, some  consideration  must  be  provided  for  c when  n > 1.  It 

n 

is  assumed  here  that  |e^|  ^ n > 1.  This  allows  the  velocity 

and  temperature  distributions  to  be  expressed  to  first  order  as 


f(n,x,t)  = fQ(ri,x)  + ejf|(n,x)  , 

e(n,x,t)  = OgCn.x)  + c^0j(ri,x)  . 


(12) 


After  substituting  Fq . (12)  into  Eqs.  (5)  and  (6),  one  obtains 

the  set  of  differential  equations  representing  the  zero  and  first- 

•> 

order  approximations  to  a nearly  quasi-steady  flow: 

Zero  Order: 


<'‘ofonn>n  + ’o‘am  * ' 'on* 


(13) 


= 2x(f,,  f,,  - f,,  f„  ) 

On  onx  onn  Ox 


onn 

Pr 


+(l+M(x))fo9„„-2.(r„^8o^ 


»0r,f0x> 


(I-!.; 


Note  that  the  viscosity  ratio,  N,  is  also  expanded  in  tlie  form 


N = No  + t:jN^, 


I 

I 

I 

I 
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•I 
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First  Order: 


(N 


o‘inn>-i  * {'  * "‘-O  ‘ohm  ' ‘ * "<■■>)  'm'ln  * (’  ' 'onn'i 


^*^^on^inx  ^in^Onx  “ ^Onn^ix  " ^inn^Ox^ 


(15) 


^ 2 ^onn  " " ^*^i^Onn^n  ’ 


0 

-if  * (l  + K(x))  - 2 (l  - M(x)) 


•'i«on  (i  - “<->)  * "‘"onSx  'mV  - %n'ix  " ®inV>  <i‘> 


+ n0 


Oti 


The  boundary  conditions  are 


^ * '>••  'o  ■ 'on  • 'i  ■ 'in  • *0  ■ '>  ®i  ■ 0 • 


n * »=  * 1,  - 0,  6„  - 0;  e,  - 0 . 


(17) 


In  a later  section  a special  version  of  these  equations  will  be  special- 
ized to  wedge  flows  where  ~ x™. 


BOUNDARY-LAYER  STABILITY 

At  present  there  is  no  theory  to  predict  the  hydrodynamic  stability 
of  an  unsteady  flow  when  unsteady  effects  are  large.  For  cases  where 
the  time-dependent  component  of  the  freestream  velocity  varies  slowly, 
one  could  develop  a theory  involving  multiple  time  scales,  but  this 
approach  has  not  yet  been  carried  out  for  the  spatial  stability  problem. 
In  this  subsection,  characteristic  times  for  the  Orr-Soimnerfeld  equation 
will  be  derived  and  used  to  establish  a heuristic  criterion  to  define 
the  region  where  unsteadiness  modifies  the  stability  analysis. 


X. 
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using this  criterion,  it  will  be  shown  that  the  stability  of  an 
unsteady  boundary  layer  can  be  studied  by  using  an  approach  that  employs 
unsteady  boundary-layer  flow  theory  and  quasi-steady  stability  analysis. 
This  approximation,  together  with  an  approximate  solution  of  the  Orr- 
Sommerfeld  equation,  will  be  employed  to  discuss  the  stability  of  the 
flows  presented  in  the  previous  section. 

Referring  to  Ref.  7 for  a heuristic  discussion  of  the  Orr-Summer feld 
equation  for  an  arbitrary  time-dependent  mean  flow,  one  will  find  tliat 
I the  constant  property  two-dimensional  version  of  the  equation  is 

h ~ ^ ~ 

(18) 

■ (is)  l♦"'’  - a‘*l  , 

where  the  superscript  prime  denotes  a derivative  with  respect  to  (y/6). 
All  velocities  are  scaled  with  R,  the  Reynolds  number,  is  based 
on  boundary-layer  thickness;  and  a is  the  wave  number. 

The  approach  that  will  be  used  in  this  subsection  will  be  to 
Investigate  the  Orr-Sommerf eld  equation  in  three  regimes:  the  outer 
Inviscld  region,  the  inner  viscous  region  near  the  wall,  and  the  crit- 
ical layer.  These  are  the  important  regions  for  a uniformly  valid 
solution  (Ref.  8),  when  matched  asymptotic  expansion  methods  are  applied. 

The  Invlscid  Stability  Equations 

In  the  limit  of  large  otR,  the  equation  reduces  to 

(<))''  - a^(|))  + (u  - cXr  - a\)  - u"(()  sr  0 . (19) 


To  determine  the  relative  effect  of  an  unsteady  freestream  veloc- 
ity, the  magnitude  of  the  first  two  terms  must  be  compared.  The  magni- 
tude of  the  ratio  of  the  unsteady  term  to  a convection  term  is  at  most 
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(6/U^ia)/0/3t)  0" 

(U  - C)  (j)" 


(20) 


where  the  magnitude  of  the  convection  term  is  approximated  by 

(u  - c)(<j)"  - a^(J))  = 0(1)  . (21) 

For  the  purpose  of  comparison  with  the  boundary-layer  quasi- 
steady parameter,  it  is  convenient  to  reduce  Eq.  (20)  to  a multiple 
of  the  ratio  of  diffusion  time  to  flow  time.  The  new  parameter  is 


a 


3U, 


ot  2 

3tU^ 


aRv  3tU, 


]_ 

aR 


(22) 


and  for  quasi-steady  flow  in  the  inviscid  stability  equation  this 
parameter  must  be  small.  Note  that  this  restriction  is  less  stringent 


than  the  boundary-layer  condition  since  aR 
boundary  layer . 


10  for  a flat-plate 


The  Wall  Viscous  Region 

Following  Ref.  8,  a scaling  appropriate  for  the  viscous  region 
near  the  wall  is 


n = y(Ra) 


F(n.t) 

(Ra) 


= <J)(y.t)  . 


(23) 


Employing  the  new  dependent  and  independent  variables  to  evaluate  the 
unsteady  term  and  the  first  viscous  term  which  is  the  dominant  viscous 
term  of  Eq.  (18),  it  can  be  shown  that  since 


6 ^ 2, , 6 3/ 

U,a  3t  - a 4))  = „ | f- 


Uj^a  3t  y nn 


2 ' 
g F 

Ra 


and 


4,"" 

aR 


nnnn 


(24) 
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If  one  estimates  that  0OF^^/U^3t)  >•  0O/3t)  (Tt  , it  can  be  con- 
cluded that  the  reduced  viscous  term  is  of  unit  order  and  the  relative 
magnitude  of  the  unsteady  term  Is  at  most 


u|a  ■ 'f  ■ O' 


Thus  the  unsteady  contribution  for  mild  unsteadiness  is  the  same  order 
of  magnitude  In  the  wall  viscous  region  as  it  Is  in  the  outer  Invlscld 
region. 

The  Region  around  the  Critical  Layer 

The  next  region  to  investigate  is  the  region  near  the  critical 
layer.  The  most  convenient  way  to  analyze  this  region  Is  to  use  the 
following  change  In  dependent  and  Independent  variables: 


n * (y  - yj.)(aR) 


♦ (y.t)  - G 


These  definitions  scale  the  unsteady  term  and  the  first  viscous 
term.  The  results  are 


.2/3  ’ 


aR  nnnn 


The  resulting  order  of  magnitude  of  the  ratio  of  these  two  terms  is 


(28) 
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T 


This  paramoter  must  be  small  for  tlie  quasi-steady  approximation 
to  be  ai'i'urate  in  the  critical  layer.  From  tills  sequence  of  estinuites 
it  is  seen  that  wlien  fj/lj-  small,  then  the  unsteady  laminar  velticlty 
profile  can  be  obtained  as  a correction  to  the  c|uasi-steady  flow 
(nearly  quasi-steady),  and  the  stability  of  the  resulting  profiles  may 
be  studied  using  classical  Orr-Sommerfeld  analysis. 

Simplified  Stability  Analysis  of  Unsteady  Flow 
over  Heated  ^ur_fac^es 

Ultimately,  a complete  analysis  of  the  stability  characteristics 
of  specific  unsteady  flows,  via  numerical  solution  of  the  Orr-Sommerfeld 
equation,  would  be  useful.  For  this  preliminary  study,  stability  analy- 
sis is  restricted  to  the  estimation  of  minimum  critical  Reynolds  number. 
Lin's  original  simplified  analysis  (Ref.  9),  and  Dunn's  subsequent  mod- 
ification for  compressible  flow  (Ref.  10)  have  been  widely  used  in 
simplified  stability  analyses.  Recently,  it  was  shown  that  a mc'dified 
version  of  this  theory  and  formulae  is  appropriate  for  water  boundary 
layers  with  variable  viscosity  (Ref.  11).  In  addition,  certain  numeri- 
cal constants  in  the  original  Dunn-I.in  formulae  have  been  revised  on  the 
basis  of  accurate  numerical  results  which  were  not  yet  available  to  Dunn 
and  Lin.  While  the  original  Dunn-Lin  theory  is  now  known  to  be  inappro- 

i 

priate  for  compressible  flow,  its  treatment  of  property  variations  is 
satisfactory.  For  completeness,  the  modified  Dunn-Lin  relations  of 
Ref.  11  are 


V 


where 


(29a) 


TTu"  U ' C 


c w 


c 


v(c) 


(29b) 
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u / n 

u'  ("v  /v  1*^  / , /*  ^ / 

w w e I A / /c  - u 
cLcJw/  \/  v/v 
\ 0 V e 


dn 


- 1 « .4 


Re  sa  28 
crit 


u’  / V \ 
— I — \ 

4 I V I • 
c \ e / 


(30) 


(31) 


The  equations,  together  with  the  boundary-layer  velocity  pro- 
i,  will  be  used  to  recalculate  I 
minimum  point  of  neutral  stability. 


files,  will  be  used  to  recalculate  Re  . « U,6,/v  evaluated  at  the 

crit  11  e 
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III.  UNSTEADY  BOUNDARY-LAYER  FXOW  OVER  HEATED  WEDGE  FLOWS 


The  approach  developed  and  discussed  in  Section  II  has  been  used 
to  investigate  the  stability  of  heated  wedge  flows  under  conditions 
of  positive  or  negative  acceleration.  This  section  describes  the 
application  of  the  general  method  and  numerical  results  to  this  class 
of  flows.  In  this  way,  the  combined  effects  of  heating,  pressure 
gradient,  and  unsteadiness  can  be  studied. 


UNSTEADY  BOUNDARY-LAYER  FLOW 

When  the  flow  is  steady  and  the  external  velocity  is  represented 
by  U^  ~ x™  (wedge  flow),  the  boundary-layer  flow  is  self-similar.  The 
boundary-layer  equations  become  ordinary  differential  equations,  and 
the  variable  M(x)  in  Eq.  (5)  is  replaced  by  the  constant  parameter  6, 
that  is  equal  to  2m/(l  + m) . For  nearly  quasi-steady  flow  and  U^  ~ x™, 
the  departure  from  a self-similar  flow  is  governed  by  the  departure  from 
a quasi-steady  flow,  and  this  is  controlled  by  the  magnitude  of  the 
parameter  Thus,  for  the  special  case  of  wedge  flows,  Eq.  (12)  re- 

duces to 


f(n,x,t)  = (fo^n)  + + m ) . 

0(n,x,t)  = ©^(n)  + e^e^(n)  , 


(32) 


where  n = qA  + m.  The  new  scaling  by  the  constant  /I  + m introduced 
into  Eq.  (32)  is  to  facilitate  comparison  to  related  works,  i.e..  Refs. 
2 and  12.  If  Eq.  (32)  is  substituted  into  Eqs.  (13)  through  (16),  the 
following  sets  of  ordinary  differential  equations  are  obtained: 

Momentum: 


* vs  * - «o>^'  - 0 • 

Vi  ■ “o'i  <”> 

- (2  - 8)  [f ■ + n ^ - i]  + (oN^e J") ' . 
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Energy : 


Pr  ^ ^0®0  ~ 


0'i'  _ _ _ no’ 

^ + fy0[  - 2(1  - B)f^0j  = - (3  - 2B)0^f J + (2  - 8)  ~ . 


(34) 


The  boundary  conditions  are  the  same  as  those  shown  in  Eq . (17). 
The  first-order  equations  for  f^  and  have  liomogeneous  boundary  con- 
ditions and  would  have  a trivial  solution  if  the  right-hand  sides  of 
Eqs.  (33)  and  (34)  were  zero.  From  this  it  is  Implied  that  the  effects 
of  unsteadiness  decrease  as  the  pressure  gradient  becomes  more  favor- 
able because  the  magnitude  of  the  inhomogeneous  term  decreases.  The 
impact  of  heating  is  not  clear. 

Equations  (33)  and  (34)  are  solved  numerically,  and  the  resulting 
parameters  are  given  in  Table  1 for  a heated  water  boundary  layer  witli 
an  ambient  temperature  of  520°R.  Moreover,  the  skin  friction  coef- 
ficient, Cj,  the  Nusselt  number,  Nu,  the  displacement  thickness,  6^, 
and  the  momentum  thickness,  6^  are  written  as 


/Red  + m)  . _ 2 - 

' (f  - f'^)  dn 


'^2,0  ^*^2,1  ’ 


m)  ^ ^ ^l'l(«>  ’ 


Nu 


■VS? 


m) 


0^(0)  + EjO|(())  ; 


(35) 


(36) 


(37) 


(3H) 


I 

I 


I 
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RKsm;i;s 

In  view  of  the  previous  disetission  and  that  of  Kef.  2,  the  restjlts 

presented  in  Table  I are  not  surprising.  In  fact,  tlie  results  for 

B = 0 and  T - T = 0 are  identical  to  those  given  in  Rei . 2.  The  new 
w e 

results  are  those  for  3 0 and  T - T ^ 0.  In  general,  we  can  con- 

w e r.  . 

i-lude  that  accelerating  flows  demonstrate  increased  skin  friction  co- 
efficient and  decreased  Nusselt  number.  These  conclusions  are  consistent 
with  those  presented  in  Ref.  2.  Further,  acce 1 erat it>n  results  in  a 
decreased  displacement  thickness  and  a small  increase  in  momentum  thick- 
ness. The  physical  explanation  for  this  is  that  acceleration  causes  an 
increase  in  the  convection  of  momentum,  and  the  wall  shear  must  increase 
to  dissipate  this  excess  momentum. 

Similarly,  for  a nearly  quasi-steady  accelerating  flow,  the  unsteady 
contribution  to  the  heat  transfer  is  affected  principally  by  tbe  thermal 
inertia,  which  resists  heat  transfer  and  leads  to  a decrease  in  Nusselt 
number.  Similar  findings  have  been  reported  in  Ref.  13.  There  is  a 
diminution  of  these  trends  as  the  pressure  gradient  bei'omes  more  favor- 
able. Naturally,  the  contrary  is  true  for  decelerating  flows. 

A word  of  caution  regarding  flows  with  adverse  pressure  gradient: 
Note  the  relatively  large  magnitude  of  fj(0)  in  comparison  to  that  of 
a flat  plate.  This  implies  that  the  unsteadiness  correction  is  larger, 
and  the  region  of  validity  of  the  two-term  expansion  may  be  more  limited 
for  adverse  gradients. 

The  stability  results  were  obtained  by  applying  Eqs.  (29)  through 

(31)  to  determine  the  critical  Reynolds  number  (based  on  displacement 

thickness,  Re  . )•  The  effect  of  flow  unsteadiness  in  the  range 
cr  it 

-.05  < e < .1  was  considered  for  a variety  of  pressure  gradients  and 
surface  overheat.  These  results  are  summarized  in  Table  2 and  are 
illustrated  in  Figs.  1 through  6. 

Using  3s  an  indication  of  stability,  it  can  be  seen  that 

for  3 = 0,  and  ~ = 0,  the  effects  of  unsteadiness  are  most  signif- 

icant. This  parameter  is  reduced  by  35  percent  when  e is  -.05,  in- 
creased by  200  percent  when  e is  +.05,  and  increases  by  almost  600  per- 
cent when  e doubles  from  ,05  to  1, 
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Table  1 

BOUNDARY-LAYER  CHARACTERISTICS 


B 

Tw  - Te 

fo<°> 

e'(0) 

q(0) 

*1.0 

-*1.1 

*2,0 

-*2.2 

0 

0 

0.46947 

1.1999 

1.218 

.7273 

.46917 

.0060 

30 

0.65699 

-1.0161 

1.4275 

.31512 

1.070 

.6232 

.4459 

.0338 

40 

0.72072 

-1.0394 

1.4919 

.31156 

1.028 

.59127 

.4378 

.0406 

60 

0.84796 

-1.0326 

1.6063 

.3094  3 

0.9525 

.53284 

.4216 

.0513 

0.1 

0 

0.58694 

0.97274 

1.080 

.44490 

.4354 

.00302 

30 

0.8062 

-1.072 

1.1533 

.45502 

0.9491 

.39775 

.4111 

.0273 

40 

0.8797 

-1.0948 

1.2046 

.44761 

.9115 

.381542 

.4029 

.0331 

60 

1.0252 

-1.1371 

1.2959 

.43780 

.8441 

.35034 

.3866 

.0423 

0.2 

0 

0.6866 

0.8208 

.9838 

.30339 

.4082 

.0014 

30 

0.932558 

-1.11379 

0.96722 

.52867 

.8632 

.28058 

.3836 

.0230 

40 

1.0143 

-1.1365 

1.0087 

.5200 

j .8287 

.2714 

.3753 

.0281 

60 

1.1752 

-1.178 

1.0825 

.507496 

1 .7668 

.2527 

.3591 

.0362 

0.5 

0 

.92776 

.57146 

1 

.8048 

.1277 

.3502 

.0127 

30 

1.2366 

-1.2002 

.62156 

.62266 

.7022 

.12988 

.3255 

.0150 

40 

1.3378 

-1.2222 

.64356 

.61370 

.6732 

. 12828 

.3174 

.0199 

60 

1.53511 

-1.2630 

.68228 

.59963 

.6215 

.1234 

.3016 

.0261 

-0.05 

0 

.400123 

1.3729 

1.313 

1.0028 

.4903 

.0077 

Table  2 

CRITICAL  REYNOLDS  NUMBER  AND  SHAPE  FACTOR  AS  A FUNCTION  OF 
UNSTEADINESS  PARAMETER  e,  SURFACE  OVERHEAT,  T - T , 

AND  PRESSURE  GRADIENT  PARAMETER  S ''  ® 


B 

Tw  - Tc 

e - 0 

c = 0.05 

e - -0.05 

c - 0.1 

e - 0 

£ - 0.05 

c - -0.05 

e » 0.1 

0 

0 

514 

1040 

332 

2940 

2.59 

2.51 

2.67 

2.44 

30 

4770 

10240 

2050 

18100 

2.40 

2.34 

.246 

2.28 

40 

7270 

13680 

3450 

21800 

2.35 

2.29 

2.41 

2.23 

60 

11100 

17900 

6190 

25700 

2.26 

2.21 

2.31 

2.16 

0.1 

0 

1380 

2750 

794 

5490 

2.48 

2.43 

2.53 

2.38 

30 

9460 

14490 

5690 

20300 

2.31 

2.27 

2.36 

2.23 

40 

12270 

17520 

7980 

23200 

2.26 

2.22 

2.30 

2.19 

60 

15570 

20470 

11170 

25710 

2.18 

2.15 

2.22 

2.12 

0.2 

0 

2880 

4800 

1750 

7610 

2.41 

2.37 

2.45 

2.34 

30 

13100 

17200 

9540 

21500 

2.25 

2.22 

2.28 

2.19 

40 

15700 

19700 

12100 

23800 

2.21 

2.18 

2.24 

2.15 

60 

18300 

21700 

14800 

25200 

2.14 

2.11 

2.16 

2.09 

0.5 

0 

7860 

9670 

6300 

11700 

2.30 

2.28 

2.32 

2.25 

30 

18900 

20800 

16900 

22700 

2.16 

2.14 

2.17 

2.13 

40 

20600 

22300 

18800 

24000 

2.12 

2.11 

2.14 

2.09 

60 

21400 

22800 

19900 

24100 

2.06 

2.05 

2.07 

2.04 

-0.05 

0 

311 

555 

216 

1681 

2.68 

2.58 

2.78 

2.48 
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From  Figs.  1 through  5,  we  observe  that  for  those  situations  in 
which  pressure  gradient  and  surface  overheat  is  small,  even  mild  un- 
steadiness (|e|  S .05)  can  have  powerful  effects  on  the  value  of  the 
critical  Reynolds  number  and  the  location  of  the  point  of  neutral 
stability.  However,  larger  values  of  favorable  pressure  gradient  and 
surface  overheat  which  are  already  characterized  by  increased  stability, 
exhibit  less  sensitivity  to  unsteadiness.  From  the  viewpoint  of  the 
boundary-layer  wall  compatibility  condition  alone,  this  result  is  sur- 
prising. However,  further  reflection  and  reconsideration  of  the  known 
effects  of  combined  heating  and  pressure  gradient  on  stability  (Ref. 

14)  suggest  that  the  situation  may  not  be  very  different  here.  For 
combined  heating  and  pressure  gradient,  it  is  known  that  the  relative 
effect  of  additional  surface  overheat  on  stability  decreases  with  in- 
creasing favorable  pressure  gradient  or  surface  overheat.  In  fact,  it 
has  been  shown  that  large  values  of  overheat  can  decrease  the  stability 
of  flat-plate  and  adverse  pressure  gradient  flows. 

Note  that  the  shape  of  the  curves  changes  when  AT  increases.  As 
suggested  above,  the  larger  the  value  of  AT,  the  less  sensitive  Re 

cr  It 

is  to  further  heating.  At  6 = .5,  for  example.  Re  . exhibits  only 

cr  it 

weak  dependence  on  further  surface  overheat  and  unsteadiness  when 
I AT  > 30°. 

Figure  5 shows  an  adverse  pressure  gradient  case  (6  = -.05)  in 
which  the  effects  of  unsteadiness  are  again  significant.  This  arises 
because  unsteadiness  dominates  the  wall  region,  and  stability  is  pri- 
marily controlled  by  details  of  the  velocity  profile  between  the  crit- 
ical layer  and  wall. 

Figure  6 shows  a collection  of  the  results  for  the  various  pressure 
gradients,  wall  heating,  and  degree  of  unsteadiness,  and  a function  of 
the  shape  factor  H = A fair  degree  of  correlation  is  demonstrated, 

indicating  that  H is,  to  a rough  approximation,  a universal  parameter  that 
can  be  used  to  estimate  the  combined  effects  of  these  several  parameters 
on  stability.  The  line  on  this  figure  is  the  result  of  exact  numerical 
computations  (Ref.  15)  for  steady  Isothermal  wedge  flows.  Similar  com- 
putations for  steady  heated  wedge  flows  in  water  (Ref.  14)  have  already 
demonstrated  that  H is  a useful  correlation  parameter  for  the  combined 
effects  of  pressure  gradient  and  heating  under  conditions  of  moderate 
heating. 


\ 


Isothermal  wedge  flows 
(Ref.  Wozzon-Okamura- 
Smith  1968) 


X P=0;  AT  = 0 
A p = 0;  AT  = 60‘’F 
□ P=0.5AT=0 
O p=  0.5AT=60‘>F 


Subscript  l-<=  0.10 

2- <=  0.05 

3- «=  0.00 

4- «=  -0.05 
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IV.  CONCLUSION 

Under  conditions  that  are  likely  to  occur  in  practice,  the  present 
nearly  quasi-steady  approach  to  the  calculation  of  unsteady  laminar 
boundary-layer  velocity  profiles  should  be  highly  accurate.  Our  analy- 
sis of  the  stability  of  unsteady  flow  is  not  yet  as  precise  as  it  will 
be  after  deeper  study  involving  multiple  time  scales.  However,  the 
heuristic  arguments  presented  here  suggest  that  there  is  an  important 
practical  flow  regime  where  classical  Orr-Sommerf eld  analysis  can  still 
be  used  to  define  the  stability  characteristics  of  nearly  quasi-steady 
velocity  profiles.  Our  analysis  and  computation,  performed  on  this 
basis,  infer  that  even  mild  unsteadiness  can  have  a powerful  impact  on 
flow  stability.  This  impact  is  particularly  strong  at  realistic  levels 
of  pressure  gradient  and  heating  in  water. 
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